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p ■ Abstract 

pH I Exotic dark matter and dark energy together seem to dominate 

• in the Universe. Supersymmetry naturally provides a candidate for 

• rH , the dark matter constituents via the lightest supersymmetric particle 

r> I (LSP). The most important process for directly detecting dark matter 

' is the LSP-nucleus elastic scattering by measuring the energy of the re- 

coiling nucleus. In the present work we explore a novel process, which 
has definite experimental advantages, that is the detection of the dark 
matter constituents by observing the low energy ionization electrons. 
These electrons, which are produced during the LSP-nucleus collision, 
may be observed separately or in coincidence with the recoiling nuclei. 
We develop the formalism and apply it in calculating the ratio of the 
ionization rate to the nuclear recoil rate in a variety of atoms includ- 
ing ^'^Ne, ^"Ar, ''^Ge, '^^Kr and ""^^^Xe employing realistic Hartree-Fock 
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electron wave functions. The obtained ratios are essentially indepen- 
dent of all parameters of super symmetry except the neutralino mass, 
but they crucially depend on the electron energy cut off. These ratios 
per electron tend to increase with the atomic number and can be as 
high as 10%. Based on our results it is both interesting and realis- 
tic to detect the LSP by measuring the ionization electrons following 
the-LSP nuclear collisions. 

Keywords: LSP, Direct neutralino search, ionization electrons, electron de- 
tection, dark matter, WIMP. 
PACS : 95.35+d, 12.60.Jv. 



1 Introduction 

From the rotational curves of objects outside the luminous galaxies it has 
been known for some time that there must be dark (non luminous) matter 
in the Universe. Furthermore the combined MAXIMA-1 [l], BOOMERANG 
121, DASI and COBE/DMR gj Cosmic Microwave Background (CMB) 
observations as well as the high precision recent WMAP data imply 
that the Universe is flat [7] and it is dominated by dark matter and dark 
energy (cosmological constant). Crudely speaking one has: 

rib = 0.05, ncDM = 0.30, i^a = o.es 

for the baryonic, dark matter and dark energy fractions, respectively. It is 
also known that the non exotic component cannot exceed 40% of the CDM 
[H]. So there is plenty of room for the exotic WIMPs (Weakly Interacting 
Massive Particles). 

Many experiments are currently under way aiming at the direct detection of 
WIMPs. In fact the DAMA experiment has claimed the observation of 
such events, which with better statistics have subsequently been interpreted 
as a modulation signal jTU]- These data, however, if interpreted as due to 
the scalar interaction, are not consistent with other recent experiments, see 
e.g. EDELWEISS dU and CDMS [12]. 

Supersymmetry naturally provides candidates for the dark matter con- 
stituents JH] 5^11- EI- ■'■^ most favored scenario of supersymmetry the 
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lightest supersymmetric particle (LSP) can be simply described as a Majo- 
rana fermion, a linear combination of the neutral components of the gauginos 
and higgsinos [T^.[ni-pi]. 

The event rates, however, are expected to be quite low and the nuclear 
recoil energies are extremely small. Thus one has to try to reduce the back- 
ground to the lowest possible level and to understand possible origins of 
backgrounds to be corrected for [23]. Anyway one has to search for charac- 
teristic signatures associated with the LSP. Such are: 

• The time dependence of the event rates, which is due to the motion of 
the Earth (modulation effect). Unfortunately, however, the modulation 
effect is small, < 2%, and for heavy targets the relevant sign depends 
on the LSP mass. So one cannot accurately predict even when the 
maximum occurs. Furthermore one should keep in mind that there are 
seasonal effects on the background as well. 

• The correlation of the observed rates of directionally sensitive experi- 
ments with the motion of the sun [221 123 12H| • In this case one has large 
asymmetries and larger relative modulation amplitudes, with charac- 
teristic directional dependence. These are, indeed, very interesting sig- 
natures, but the expected event rates are suppressed by a factor of at 
least Att relative to the standard rates. Furthermore, at present, the 
experiments are able to measure the direction of nuclear recoil, but it 
will be extremely difficult to determine the sense of direction j2H]- 

• Inelastic excitations to low lying excited nuclear states have also been 
considered [301 1^ . In this case one can overcome the difficulties en- 
countered in the recoil experiments (energy threshold, quenching fac- 
tors etc) by detecting the de-excitation 7-rays. The predicted branching 
ratios, due to the spin induced DM cross sections, in favorable circum- 
stances can be as high as 10%, which is a very encouraging result. 

• Detecting the dark matter constituents by observing the low energy 
electrons, which follow the ionization of the atom during the LSP- 
nucleus collision fIE\ . 

The last possibility may be realized with the technology of gaseous TPC de- 
tectors In fact the WIMP-nucleus scattering leads: (i) to nuclear 
recoil without atomic excitations and (ii) to nuclear recoil with atomic excita- 
tions. So far most CDM searches have been made by the inclusive processes. 
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(i) and (ii), employing solid detectors. We propose that the produced elec- 
trons in atomic excitation (ii) should be studied by exclusive measurements. 
In a previous calculation we have estimated the ratio of the ionization 
rate (per electron) by considering the simple target ^°Ne assuming hydro- 
genic electron wave functions. 

In the present paper we extend these calculations by considering a number 
of experimentally interesting targets. To this end realistic Hartree-Fock elec- 
tron wave functions have been employed. The observed ratio of the electron 
rate (ii) for one electron per atom, to be denoted below as per electron, di- 
vided by the standard nuclear recoil rate (i), depends, of course, on the low 
energy cut off for both electrons (tenths of keV) and nuclear recoils (10-20 
keV) as well as the quenching factor for recoils. They can reach It is found to 
increase as the atomic number and it can get as high as 10% for a reasonable 
choice of these parameters. The experimentally interesting quantity the ra- 
tio per atom, meaning that one considers all electrons in the atom, obtained 
from the previous one by multiplying it with Z . This ratio for medium and 
heavy nuclei may even become larger than unity. Thus this method can be 
widely used for any target nuclei with adequate electron detection capability. 
The present paper is organized as follows: The theoretical framework is pre- 
sented in Sec. El the main results are shown in Sec. El and a brief summary 
is given in Sec. HI 

2 The essential ingredients of our calculation 

The electron ionization due to the LSP collisions can take place via two 
mechanisms: 

1. The electrons are ejected during the LSP-nucleus collisions via the stan- 
dard EM interactions. 

2. The electrons are ejected due to the LSP-electron interaction. This 
interaction is obtained in a fashion analogous to the LSP-quark inter- 
action. In the case of the Z and Higgs exchange one replaces the rel- 
evant quark couplings (23 by the electron couplings. In the s-fermion 
mediated mechanism the s-quarks are replaced by s-electrons. In fact 
this mechanism has recently been proposed as providing a possibility 
of LSP detection in electron accelerators (3^]. In the present case, how- 
ever, one can get reasonably high electron energies, of a few hundreds 
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of electron volts, only if one exploits electrons bound in nuclei. Then 
the kinematics is similar to that of the previous mechanism. 

Since the branching ratios of the second case depend on the SUSY parameter 
space and one does not enjoy the effect of coherence arising from the scalar 
interaction, we will concentrate here in the first mechanism. 

2.1 The formalism 

The differential cross section for the LSP nucleus scattering leading to the 
emission of electrons in the case of non-relativistic neutralino takes the form 



where v, and are the oncoming LSP velocity, energy and momentum 
distribution, while q is the momentum transferred to the nucleus. M is the 
invariant amplitude, known from the standard neutralino nucleus cross sec- 
tion, T and k are the kinetic energy and the momentum of the outgoing 
electron and e^; is the energy of the initial electron, which is, of course, neg- 
ative. 0„im(k) is the Fourier transform of the bound electron wave function, 
i.e., its wave function in momentum space, pni is the probability of finding 
the electron in the n, I orbit. In the expression above and in what follows 
our normalization will consist of one electron per atom, to be compared with 
the cross section per nucleus of the standard experiments. 

In order to avoid any complications arising from questions regarding the 
allowed SUSY parameter space, we will present our results normalized to the 
standard neutralino-nucleus cross section. The obtained branching ratios 
are essentially independent of all parameters of supersymmetry except the 
neutralino mass. 

With these ingredients we find that the ratio of the cross section with ion- 
ization divided by that of the standard neutralino-nucleus elastic scattering, 
nuclear recoil experiments (nrec), takes the form 
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'^n^Px^ = 2^r'^S, is the momentum q transferred to the nucleus and F{q) 
is the nuclear form factor. For a given LSP energy, the outgoing electron 
energy lies in the range given by: 



< T < -^T^ + eni. 



(5) 



Since the momentum of the outgoing electron is much smaller than the mo- 
mentum of the oncoming neutralino, i.e., K 1^ the integration over can 
be trivially performed. Furthermore, if the effect of the nuclear form factor 
can be neglected, the integration over ^ can be performed analytically. Thus 

we get 
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Otherwise the angular integrations can only be done numerically. Finally 
one must integrate numerically the above expression over the electron spec- 
trum to obtain the total cross section as a function of the neutralino energy 
(velocity). The folding with the LSP velocity will be done at the level of the 
event rates (see sec. 12. 4j) 
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2.2 The role of the nuclear form factor 

The nuclear form factor F{u) entering in Eq. Q has the general form 

Z N 

Fiu) = jFziu) + -Fr^iu) (7) 

The proton and neutron form factors (-Fz(m) and F]\f{u) respectively) are 
calculated by the Fourier transform of the proton density distribution pz{r) 
and neutron density distribution pAr(r) respectively. In the present work 
p^(r) and PAr(r) are constructed using harmonic oscillator wave functions. 
More specifically the densities are defined as 

pz,^(r) = 5:0*(r)0,(r) (8) 

i 

where the sum run all over the proton ( neutron) orbitals 0i(r). In general 
some of the nuclei considered in this work may be deformed. Thus, strictly 
speaking, the densities depend not only on the radial coordinate r but on the 
angle 6 as well. In the present work, however, we are mainly interested in 
branching ratios, which depend very little on the nuclear form factor. Thus 
to a good approximation the nuclei can be treated as spherical. Thus in the 
case of the heaviest nuclear system considered in this work, namely ^^^Xe, 
for a threshold energy of Eth = 0.2 keV the effect of the form factor is a 
reduction of the branching ratio by less than 7%. For this reason we have 
included only one stable isotope, namely the most abundant one, which is 
indicated on the figures. 

In view of the above the form factor is given by the relation: 

FzMq) = I e'^'pz,N{r)dv (9) 

The nuclear form factor of ^°Ne is taken from Ref. |SZI, where the nuclear 
wave functions were obtained by shell model calculations using the Wilden- 
thal interaction. The nuclear form factors of the heavier nuclei were for 
simplicity calculated by using spherical harmonic oscillator wave functions. 
This is a good approximation since, as it has already been mentioned, the 
ratio of the two processes we is insensitive to the details of the nuclear wave 
functions. 
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2.3 The bound electron wave functions 



In the present work we consider very accurate spin-independent atomic wave 
functions obtained by Bunge et al ^j, by applying the Roothaan-Hartree- 
Fock method (RHF) to calculate analytical self-consistent-field atomic wave 
function. In principle one should employ the relativistic electron wave func- 
tions, especially for the inner shell electrons of large Z atoms. Anyway for 
our purposes the above wave functions are adequate, since the inner shell 
electrons contribute only a small fraction to the rate. In this approach the 
radial atomic orbitals Rne are expanded as a finite superposition of primitive 
radial functions 

Rniir)=J2C,„iS,iir) (10) 
j 

where the normalized primitive basis Sje{r) is taken as a Slater-type orbital 
set, 

Sji{r) = N,,r'^^'-^e-^^'' (11) 
where the normalization factor Nj^ is given by 

N,, = (2Z,,)(«^-+i/2)/[(2n,,)!]i/2 (12) 

and rtjf^ is the principal quantum number, Zji is the orbital exponent, and £ 
is the azimuthal quantum number. 

The atomic wave functions in the momentum space are related to the 
coordinate wave functions by 

0n£m(k) = j e-'^^'ct^nhnWr (13) 

where nlm denote the usual quantum numbers characterizing atomic states. 
The radial momentum-wave function defined as 

^nimiX) = i-tYRneik)YeU^k) (14) 
is related to the radial wave function in coordinate space through (see also 

m) 

Rni{k) = J-j^ r'Rni{r)j,{kr)dr (15) 

where ji{kr) is a spherical Bessel function. The radial wave functions in 
momentum space Rni{k) are written as 

Rni{k)=Y.CjnlS,,{k) (16) 
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in term of the RHF functions Sje{k) in momentum space, related to Sje{r) 
through 

SAk) = \ - / r^S,,{r)je{kr)dr (17) 
V vr JO 

The needed binding energies were taken from [lUj . 

2.4 Folding with the LSP velocity distribution 

The above resuhs were derived assuming a definite LSP velocity. In practice, 
however, the LSP obeys a velocity distribution. The actual velocity distri- 
bution may be complicated, i.e. asymmetric and varying from one point to 
another. In our vicinity it is commonly assumed that LSP obeys a simple 
Maxwell-Boltzmann velocity distribution with respect to the galactic center 
P5| . namely 

^<''' - (^^^"'"'""' 

with vq = 220Km/s. This velocity distribution does not go to zero at some 
finite velocity. So on this one imposes by hand an upper velocity bound 
(escape velocity), Vesc = 2.84vo- One then must transform this distribution 
to the lab frame, v ^ v + Vq, since Vq is also the velocity of the sun around 
the center of the galaxy. We will not be concerned here with the motion of 
the Earth around the sun, i.e., the modulation effect [2HI 

Folding both the numerator and denominator of Eq. Qwith the LSP 
velocity distribution, after multiplying each with the LSP flux 

p(0) m 



Anip 

we obtain the differential ratio with R^, the rate for the ionization, on 

the form 



1 dRe da{T) 



R dT a„ 



^Pnt\4>nfi\J'2,meT)\'^me^2meT 



V -V K sinh 2t;/t;o rft; , , 

jt^rna. ^^2g-i.V"o'sinh 2i;/t;o rfw 



Ivir, 

where the numerator factor and the denominator factor D are given re- 
spectively 



2 /■'^3 

N = — q,F\q,)dq,, q, = qb (20) 
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D = \ r qbF\qb)dqb, Qb = qb (21) 
af Jo 



where 



«. = -^. «= = isr«- -3 = ^(i + \/r^) (22) 

and b is the width of the harmonic oscillator potential. 

In Eq. (fTUI) Vmax is the escape velocity, while Vmin is given by Eq. (0), 

/ 2(T-6„,) 

V /^r 

while for the nuclear recoils: 



j2AmpQth 

'^^min ~ v24j 

It is clear that the fraction in Eq. (|19p is a decreasing function of the outgoing 
electron energy T as well as the binding energy — e„£ of the initial electrons. 
For T in the sub-keV region, where the differential rate becomes sizable, Vmin 
is determined essentially by the energy of the bound electron. This restriction 
on the velocity distribution is a source of suppression of the rate associated 
with the inner shell electrons in the case of heavy targets, in particular, the 
Is electrons. In fact for Xe for T = this fraction reaches the maximum 
values of 0.13,0.75,0.76 and 0.95 for Is, 2s, 2p and 3s respectively. For 
lighter systems this kinematical suppression is not important even for the Is 
orbits. Further suppression comes, of course, from the behavior of the bound 
electron wave function. 



3 Results and Discussion 

We begin our discussion by focusing on the differential rate. It is useful to 
divide it by the total nuclear recoil rate, i.e. in order to make our 

results independent of the details of the SUSY parameter space and related 
uncertainties. From previous work [23] on ^''Ne with hydrogenic wave func- 
tions we know that the differential rate sharply picks at low electron energies. 
This result is confirmed by our present calculation employing realistic elec- 
tron wave functions. This is exhibited by considering three different atoms 
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(^°Ar, '''^Ge and ^^^Xe) for the typical LSP mass = 100 GeV (see Figs. 
miSI). If we then integrate the differential rate with respect to the elec- 
tron kinetic energy T, we obtain the relevant event rate ratio The thus 
obtained results for a number of nuclei are shown in Fig. They are pre- 
sented as a function of the threshold energy for electron detection Eth for 
LSP mass = 100 GeV. We also present them as functions of the LSP 
mass m^, for Et^ = 0.2 keV (Fig. We clearly see that the results are 

very sensitive to the threshold energy. It is also clear that the heavier tar- 
gets are favored. It is encouraging that branching ratios per electron of 10% 
are possible, if one can reach threshold energies as low as 200 eV, which is 
feasible for gas targets . It should be mentioned that in obtaining these 
results we assumed threshold effects for the ionization electrons, but not for 
the nuclear recoils. The nuclear recoil rate is less sensitive to the threshold 
in the sub-keV region. 

As is well known the recoil experiments suffer from the effect of quenching. 
The idea of quenching is introduced, since, for low energy recoils, only a 
fraction of the total deposited energy goes into ionization. The ratio of the 
amount of ionization induced in the gas due to nuclear recoil to the amount 
of ionization induced by an electron of the same kinetic energy is referred to 
as a quenching factor QufaciQ)- This factor depends mainly on the detector 
material, the recoiling energy Q as well as the process considered |^T]. In 
our estimate of QufaciQ) we assumed a quenching factor of the following 
empirical form motivated by the Lidhard theory j^-|42j: 

, ri ~ 0.256 , ra ~ 0.153, (25) 

where Q is the nuclear recoil energy. The quenching factor QufadQ) is shown 
in Fig. m Since, as we have mentioned, we expect a threshold energy of 10 
keV and the quenching factor becomes ~ 0.3 above this energy, the branch- 
ing ratio will increase by another factor of three, if the quenching is taken 
into account. 

The number of electron events divided by the nuclear recoil events is obtained 
by multiplying the above branching ratios with the atomic number Z. This 
enhancement makes it possible to detect electron events not only with TPC 
detectors but nuclear double beta decay detectors ) sensitive to detection 
of low energy electrons ~ IkeV. 

With this number of events one wonders, if one can detect X rays formed 
after the inner shell electron hole, created by the LSP- nucleus collision, is 



QfaciQ) = ri 



IkeV 
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filled by one of the outer electrons, with or without coincidence measurement 
of the outgoing particles (electrons and/or recoiling nuclei). In fact we find 
that some inner shell orbits, in particular 2p and 3d, can make a sizable con- 
tribution. It thus seems that such X rays can be exploited. Discussion of this 
phenomenon as well as the background of the Auger electrons involving the 
inner shell electrons [221 will be discussed in detail in a separate publication. 
So far, the ratio of the electron rate to the recoil rate is normalized to one 
electron per atom. In practice, the total electron rate for given detector mass 
is obtained by multiplying the above rate per electron by the atomic number 
Z . Then the ratio of the total electron rate to the recoil rate is plotted as a 
function of the threshold energy in Figs El El Here in addition a quenching 
factor of ~ 0.3 is used for the nuclear recoil. The thus obtained ratio is 
very impressive, even at an electron threshold higher than 0.5 keV. Thus it 
is interesting to search for electrons with the high-sensitivity low-threshold 
detectors such as those used for double beta decays discussed in recent review 
papers [IH] 

4 Concluding Remarks 

We have explored a novel process for the direct detection of dark matter, that 
is by exclusive observation of the low energy ionization electrons produced 
during the LSP-nucleus collision. We have applied our formalism in a variety 
of atoms employing realistic electron wave functions and nuclear form factors. 
We have focused in noble gases, since we expect that low energy electrons 
can better be detected by gaseous TPC detectors. Our results, however, can 
be applied, if necessary, in other targets, e.g. those employed in neutrinoless 
double beta decay, since we find that the details of the electronic structure 
are not very important. 

To minimize the uncertainties arising from SUSY and to make the presenta- 
tion as simple as possible, we have calculated the ratio of the rate associated 
with this process divided by the rate for the standard nuclear recoils. We 
find that this ratio per electron slightly increases with the atomic number. 
It depends appreciably of the electron and the nuclear recoil energy cut offs. 
It is quite important to make the electron detection energy cut off as low as 
possible. Branching ratios (per electron) of about 10% can realistically be 
expected for reasonable choices of these threshold energies. Thus the ratio 
of the number of electron events per atom divided by the number of recoils, 



12 



obtained from the previous ratio after multiplication with the atomic number 
can even become larger than unity. Thus, since the ionization electrons 
extend to the keV region, their contribution has to be taken into account in 
the the standard nuclear recoil experiments as well. 

As it has been previously shown the background problem are no worse 
than those appearing in the standard recoil experiments. They can be further 
reduced by coincidence experiments, by measuring the ionization electrons 
in coincidence nuclear recoils. Finally it should noted that the inner shell 
electron holes, created during the LSP-nuclear collision, can be filled by emit- 
ting X rays and/or Auger electrons. Such issues, which are of experimental 
interest, will be discussed in detail elsewhere. 

In view of the above considerations it is realistic to search for WIMPs by 
measuring ionization electrons. 
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Figure 1: The differential rate for ionization electrons, divided by the total 
rate associated with the nuclear recoils, as a function of the electron energy T 
(in keV) for various atoms. The results exhibited were obtained for a typical 
LSP mass = 100 GeV by including a nuclear form factor, but without 
recoil threshold effects. 
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Figure 2: The same as in Fig. except that we have now separated the plot 
in various energy intervals, so that the energy dependence at higher energies 
becomes visible. 
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Figure 3: (a) The total ionization rate per electron divided by the standard 
nuclear recoil rate as a function of the electron threshold energy. The results 
exhibited were obtained for a typical LSP mass rriy. = 100 GeV by including 
a nuclear form factor, but no threshold effects on recoils, (b) The same 
quantity plotted as a function of the mass of the LSP. The results were 
obtained for a threshold energy Eth = 0.2 keV 
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Figure 4: The quenching factor QufaciQ) as a function of the recoil energ 
Q in keV . 
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Figure 5: The same as in Fig. except that now the ionization rate Rg 
per electron is multiphed by the atomic number Z to obtain the rate per 
atom. In other words we plot the relative ionization rate ZRe per atom with 
respect to the nuclear recoil rate, i.e. r = Z^. The quenching factor of 
about 1/3 has also been included in the denominator. More precisely we use 
(1/3) R for the denominator (standard recoil) to incorporate the reduction 
of the recoil rate with a typical threshold and quenching factor (see Fig. HJ. 
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Figure 6: The same as in Fig. except that now the electron ionization 
rate per electron has been multiplied by the atomic number Z (the relative 
event rate per atom is r = Z^). The quenching factor of about 1/3 has also 
been included in the denominator (standard nuclear recoil). 



23 



